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1 Introduction 



Non-Hermitian quantum mechanical problems have attracted much attention 
recently. The main reason for this is that the energy spectrum of a number of 
complex potentials turned out to be real (at least partly), which contradicted 
the usual expectations regarding non-Hermitian systems. Strangely enough, 
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the first examples for complex potentials with real spectra were found using 
numerical techniques [1]. This unusual behaviour of the energy spectrum was 
attributed to the so-called PT symmetry, i.e. the invariance of the Hamil- 
tonian with respect to the simultaneous space ( P) and time ( T) reflection. 
For one-dimensional potential problems this requires [V(— x)]* = V(x), which 
implies that the real component of the potential must be an even function of 
x, while the imaginary component has to be odd. After the first examples, 
further ones have been identified using semiclassical [2], numerical [3] and 
perturbative [4] methods, and a number of exactly solvable PT symmetric 
potentials have also been found, mainly as the analogues of Hermitian (real) 
potentials [5-7]. 

It was also noticed that PT symmetry is neither a necessary, nor a sufficient 
condition for having real energy spectrum in a complex potential. It is not a 
necessary condition, because there are complex non- PT symmetric potentials 
with these properties: some of these are complex supersymmetric partners of 
real potentials [8], while some others can be obtained by merely shifting PT 
symmetric potentials along the x axis, which (formally) cancels PT symme- 
try, but obviously does not influence the energy spectrum. Neither is PT 
symmetry a sufficient condition, because complex-energy solutions of such po- 
tentials are also known, and since in this case the energy eigenfunctions cease 
to be eigenfunctions of the PT operator, this scenario has been interpreted 
as the spontaneous breakdown of PT symmetry [1]. No general condition has 
been found for the breakdown of PT symmetry, but it has been observed that 
it usually characterizes strongly non- Hermitian problems [1,9-11]. 

Obviously, the lack of Hermiticity raises questions about the probabilistic 
interpretation of the wavefunctions (probability density, continuity equation), 
and in general, about the definition of the norm and the inner product of 
the eigenvectors of the non-Hermitian Hamiltonian. It has been suggested, for 
example, that the ifj 2 (x) quantity should replace [^(x)) 2 in the definition of 
the norm [12]. For unbroken PT symmetry this expression coincides with the 
ijj(x)il>*(—x) quantity used in the definition of the pseudo-norm [13], which is 
obtained from the modified inner product P\ipj). This redefinition of the 
inner product was found to lead to the orthogonality of the energy eigenstates, 
but it also resulted in an indefinite metric, replacing the usual Hilbert space 
with the Krein space [14]. Efforts have been made to restore the Hermitian 
formalism using projection techniques [13,14]. 

Exactly solvable examples can be extremely useful in the understanding of the 
unusual features of PT symmetric problems and the underlying new physical 
concepts. For example, by the continuous tuning of the potential parame- 
ters through critical values the mechanism of PT symmetry breaking can be 
studied; conditions for positive and negative values of the pseudo-norm can 
be identified, etc. Although a number of exactly solvable PT symmetric po- 
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tentials have been identified, these questions have been addressed only in very 
few cases. 

The Scarf II (sometimes also called hyperbolic Scarf or Gendenshtein) po- 
tential is in a special position among exactly solvable potentials. This shape- 
invariant potential [15] is defined on the whole x axis, it has no singularity 
at x = 0, and in contrast with most other shape-invariant potentials, it can 
be turned into a PT symmetric form without regularizing its singularity by 
means of an x — > x + ie imaginary coordinate shift [6,16]. Therefore it is not 
surprizing that it became a "guinea pig" of testing PT symmetry on a solv- 
able example. It has been associated with the sl(2,C) [17], su(l,l)~so(2,l) 
[18] and so(2,2) [19] potential algebras, and it has also been observed that 
its PT symmetric version has a second set of bound states, which appear as 
resonances in its Hermitian version [17,18]. This mechanism of doubling the 
bound states is essentially different from the one arising from the cancellation 
of singularities at x = by the imaginary coordinate shift. This potential 
is also known to have (purely) real and (purely) complex energy spectrum, 
depending on the relative strength of its real and imaginary component [9], 
and since the two domains can be connected with a continuous tuning of the 
parameters without crossing a singularity, it is a perfect example to illustrate 
the breakdown of PT symmetry. 

It would also be a suitable example to illustrate the modified definition of the 
inner product and the behaviour (sign) of the pseudo-norm, and the other 
implications for the use of the Krein space instead of the usual Hilbert space, 
however, there is a major obstacle: the evaluation of integrals involving them 
could not be calculated as yet analitically, except for the ground state [20] . In 
fact, even an explicit proof of the orthogonality of the bound states is missing, 
both for the Hermitian [21], and the PT symmetric case [9]. In the latter case 
only indirect proof has been given for orthogonality of some states [9]. 

Obviously, the evaluation of integrals containing the energy eigenfunctions 
of the Scarf II potential is essential to complete the analysis of this perfect 
example for PT symmetry, i.e. to study the behaviour of the pseudonorm 
and to prove the orthogonality of the eigenstates in general. In this Letter we 
present a method to evaluate these integrals for the first time, both for the PT 
symmetric and the Hermitian version of the Scarf II potential, and use these 
results to illustrate the mechanism of PT symmetry breaking. In particular, 
we prove the orthogonality of the eigenstates and derive the normalization 
constants (the norm) both for the PT symmetric and the Hermitian version 
of this potential, and also analyse the PT symmetric conventional 
complex potential, using the standard (Hilbert space) definition of the inner 
product. 
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2 The general form of the Scarf II potential 



Here we follow the notation of Refs. [6,22] to discuss the Scarf II potential 

V(x) = - 



cosh (x) 




2isinh(x) ( (3 + a\ { (3 — a\ . . 

+ cosh\x) \~2~) \~2~) ' { ) 

The bound-state energy eigenvalues are 

2 



(2) 



while the corresponding wavefunctions 

ij>W\x) = C^\l -isinh(a;))f + i(l + isinh(x))i + ^Pi a ' /3) (isinh(a;)) (3) 

are expressed in terms of Jacobi polynomials [23] and are normalizable if 
n < -[Re (a + (3) + l]/2 holds. 

In the Hermitian case a and (3 are complex and satisfy a* = (3: a = —s — 
\ — iA, (3 = — s — | + iA [21,24]. In this case only one regular solution exists. 
With arbitrary a and (3 obviously, the general complex version of the Scarf II 
potential is obtained. 

The Scarf II potential can be made PT symmetric if a* = ±a and f3* = ±/5 
holds [6], i.e. if a and (5 are both either real or imaginary. In order to have real 
energy eigenvalues both a and (5 have to be real, while to have complex bound 
state spectrum, i.e. in the case of spontaneous breakdown of PT symmetry 
one of them has to take an imaginary value [22]. If both a and /3 are imaginary, 
then there are no bound states. Here we assume that (5 is real, and a can be 
real or imaginary, depending on whether the PT symmetry is unbroken or 
broken. This choice does not restrict the generality of the problem, since the 
roles of a and (3 can easily be reversed (see (21) in the Appendix). 

For the Scarf II potential the breakdown of PT symmetry takes place when 
the strength of the imaginary potential component exceeds a certain limit 
depending on the strength of the real potential component, as described in 
Ref. [9]. This condition corresponds exactly to taking imaginary values for a 
instead of real ones (see e.g. Ref. [22] for the details), so a smooth transition 
over the critical point can be achieved by moving a to zero along the real axis 
and then continuing along the imaginary axis. 
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In the PT symmetric case there are two sets of normalizable solutions [6,17,18] , 
which carry the upper indexes (cx,/3) and (—a,/3) in (3). Obviously, (1) is not 
sensitive to the + or — sign of a. In the notation of Ref. [20] the two solutions 
corresponds to quasi-parity q — +1 and — 1. This sign difference results in two 
distinct energy eigenvalues in (2), which form a complex conjugate pair when 
a is imaginary, i.e. in the case of broken PT symmetry. In this case the PT 
operation transforms the two solutions into each other, while in the unbroken 
symmetry case the two solutions are eigenfunctions of the PT operator. 

In the following sections we are going to evaluate integrals containing the stan- 
dard and PT symmetric inner product of wavefunctions of the type ^ ±a,/3 ^ (x) 
(3). The technical details of the calculations can be found in the Appendix. 



3 The PT symmetric inner product and the pseudo-norm 

Let us consider the PT symmetric inner product [13,14] of two solutions of 
the type (3) 

oo 

J ^\X)[^\-X)]*dx . (4) 

— oo 

According to our choice, (3 is real and 5 can be ±a, depending on whether we 
calculate the PT symmetric inner product of states with the same or different 
quasi-parities (5 — a and 5 = —a, respectively), furthermore, a can be real or 
imaginary, depending on whether the PT symmetry is unbroken or broken. 
Applying (21) and (22) presented in the Appendix, we can write 

= c^[^n-irq^ M . (5) 

This formula together with (29) has significant implications regarding the 
PT symmetric inner product (4). First note that whenever a = —5* holds, 
the integral vanishes due to the presence of the sin[7r(a + 5*)] term in (29). 
This corresponds to either a — 8 with imaginary a, i.e. the inner product 
of wavefunctions of the same type (same quasi-parity) in the broken PT 
symmetry case, or a = —5 with real a, i.e. the inner product of two different 
type (different quasi-parity) wavefunctions in the unbroken PT symmetry 
case. So we can conclude that the two states are orthogonal in these situations. 

Now let us consider the cases when a ^ —5*. The first case is 5 = a with real 
a (unbroken PT symmetry). Substituting this S in (34) (and remembering 
that (3 is real) we get 
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l nl 



2 a+/3+2 S i n ( 7ra ) sin(7T/3) 

a + (3 + 2n + l sin[7r(a + (3)] 



x 



' a + (3 + 2n N 


-i 

1 




+ (3 + 2n^ 








n ) 



(6) 



This proves directly the orthogonality of the states of the same type (i.e. those 
with the same quasi-parity) for n ^ I when the PT symmetry is unbroken, 
and gives a closed formula for the pseudonorm for n — I. Previously this 
pseudonorm was known only for the ground state n = [20], while the or- 
thogonality of the eigenf unctions was proven only indirectly [14,9]. This latter 
proof rests on the equation 



oo 

(E n — E*) J i/> n (x)tf(-x)dx = , 



(7) 



which is the equivalent of the equation proving the real nature of the energy 
eigenvalues for Hermitan systems. In the case of unbroken PT symmetry E n 
and Ei are real and they are not equal, consequently the integral in (7) has to 
vanish. 

The only remaining case is 5 = —a with imaginary a, when 5* = a holds again. 
This case gives us the overlap of eigenstates belonging to different quasi-parity 
in the broken PT symmetry case. It turns out that the I^' 13 ^^ overlap has the 
same form as (6), except that \C^\ 2 has to be replaced with C( a> ^[C, ( " a,/3) ]*. 

Let us summarize the results for the different cases. 

• Unbroken PT symmetry (a real), same type wavefunctions: J^' /3 ' Q ' ) is di- 
agonal in n and /, as seen from Eq. (6). To extract more information, we 
can rewrite Eq. (6) in a somewhat different form, after eliminating the sine 
functions from the formulas by combining them with some gamma functions 
via r(x)r(l — x) = 7r/ shi7r:r: 

= u-ir^r )|2 — n-g-P-n) 



-a-(3-2n- l)n\ T(-a - n)T(-(3 - n) 



Due to the condition for having bound states, i.e. n < — [Re (a) + (3 + l]/2, 
if a is real, then in Eq. (8) every term is positive, except (—1)™ which 
alternates, and [r(— a — n)T(—l3 — n)] _1 , which is real, but its sign depends 
on the relative magnitude of a, (3 and n. Except for extreme values of a 
and (3 the argument of the two gamma functions is positive for the first few 
n's, so then the alternating (— l) n factor determines the sign of the pseudo- 
norm, but as n reaches —a and/or —(3, this regular pattern changes. The 
results of this case are new, except for n — 0. 
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• Unbroken PT symmetry (a real), different type wavef unctions: = 
0, due to shi7r(a — a*) = in (29). This was proven indirectly by (7) [14,9]. 

• Broken PT symmetry (a imaginary), same type wavefunctions: I^'^'^ = 0, 
due to sin7r(a — a*) = in (29). This was proven indirectly by (7) [14,9]. 

• Broken PT symmetry (a imaginary), different type wavefunctions: I^f'^' 

is diagonal in n and /, as seen from Eq. (6). But in this case it seems 
that for n — I there can be two different wavefunctions which are not 
orthogonal, in general. Equation (8) holds for this case too, except for a 
change in the term containing the normalization constants, as discussed 
before. This non-orthogonality of two different states is a new feature of 
PT symmetric problems, which in this case appears only when the PT 
symmetry is broken. This unusual result seems to be supported by Eq. 
(7): when the PT symmetry is broken, the energies of the two states with 
the same principal quantum number n but with different quasi-parity are 
complex conjugate to each other, so the zero value of (7) is secured by the 
energy term, and the integral need not be zero. 



4 The Hermitian Scarf II potential and the normalization coeffi- 
cients of the wavefunctions 



As discussed in the Introduction, the normalization coefficients of the wave- 
functions have not been determined as yet, due to the involved mathematics 
[21,9]. Here we use our new approach and show that it can also be applied to 
evaluate integrals for the conventional Scarf II potential. 

Let us now denote the (standard) inner product of the states (3) as 

oo 

= J ^\x)[4^\x)Y&x (9) 

— oo 



Remember that in this case there is only one set of bound-state eigenfunctions, 
and that in this case we have a* = (3. It is easy to show that 



Note that 8* = a* = (3 and 7* = ft* = a holds now, so the conditions under 
which the general integral can be reduced to a simple form are again satisfied. 
The final result is rather similar to (8) obtained for the PT symmetric case, 
except that the alternating (— l) n factor is now missing: 

Knl ~ 6ninlL/n 1 (-a - p - 2n - l)n\ T(-a - n)r(-(3 - n) ' W 
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It is easy to verify that this integral is always positive as it should. For this 
we recall that a* = (3, so a + (3 is real. We also recall the condition for bound 
states n < — (Re[a + /3] + l)/2, so the only terms which can introduce negative 
(and complex) quantities are the two gamma functions in the denominator. 
However, their product is now Y(—a — n)[T(— a — n)]* = a — n)\ 2 > 0, 
so we conclude that K^f^ > 0. Equation (11) thus determines the normaliza- 
tion coefficients of the bound-state wavefunctions of the conventional Scarf II 
potential for the first time: 



r(-a - n)V(-P - n)(-a -0-2n- l)n\ 
r(— a — j3 — n)n 



1/2 



(12) 



5 The PT symmetric Scarf II potential as a conventional complex 
potential 



As we have discussed in the Introduction, since the PT symmetric Scarf 
II potential is defined on the real x axis, it can also be considered as an 
ordinary complex potential, and the same overlap and normalization integrals 
can also be evaluated for it too. This means using the wavefunctions of the 
PT symmetric problem in inner products of the type (9). As expected for 
a complex potential, the states cease to be orthogonal, as it can be shown 
by direct calculation of integrals of the type (10) with real values of (3 and 
real or imaginary values of a. Technically this is reflected by the fact that the 
conditions 5* = a and 7* = j3, which were used in the Appendix to bring the 
sums to a closed form are not satisfied. However, some interesting results can 
be obtained calculating some diagonal matrix elements and overlaps. 

Starting from the one-dimensional Schrodinger equation for a generic complex 
potential 

- iP"(x) + [U(x) + iW(x)]i/>(x) = (E R + iEjMx) (13) 



for normalizable eigenfunctions ip(x), we can derive [25,26] the following rela- 
tion connecting the imaginary component of the potential and of the energy 
eigenvalue: 



This relation can be rather useful in demonstrating the PT symmetry break- 
ing mechanism: tuning the potential parameters from the domain of unbroken 
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PT symmetry to the domain of symmetry breaking should result in the ap- 
pearance of a non-zero value of Ej. For the Scarf II potential this means shift- 
ing a along the real axis up to zero, and then continuing along the imaginary 
axis. 



Specifying (14) to (1) and (3), we get 

roo ; (a,/?) ( \(P 2 -a 2 ) sinh x y ( a ,/j) , s i * i j( a fl 

" " " r oo 4 Q ' /3) (^)[4 Q ' /3) (^)]*da; ' " 



for the integrals, where we assume that ip^\ x) is a normalizable eigenfunc- 
tion of the PT symmetric Scarf II potential, i.e. f3 is real, and a is either 
real or imaginary, depending on whether the PT symmetry is unbroken or 
broken. The integral in the denominator can be expressed implicitly in a form 
similar to the diagonal version of Eq. (4): 



1 " 1 sin[vr(/3 + (a + a*)/2)] 



x 



■ n + a] [ n + fl I " , n + a*\ n + 8 

E(-!) m E(-i) m 

m=o \ m I \ n — m I m'=o \ w! \ n — m! 



r(2±^ + n + 1 - m + m')r(^ + n + 1 + m - m 1 ) 

+ B + 2n + 2) 



± y 2 r it T J- — in t ; J- v 2 r "< J- ^ — "t; , , 



This integral can also be expressed in terms of the Qi"'^' 7 ' 5 * 1 quantity (22). 
However, in this case we have 7 = a and 5 = 3, so the conditions (31) 
for reducing the implicit summed expression into a closed formula are not 
satisfied. Nevertheless, it is technically not too involved to evaluate (16) for 
the first few values of n. 

A similar expression holds also for the integral appearing in the numerator, 
which can be brought to a sum form using (28): 



j(a,/J) = 1(^2 _ a 2^ c ( a ,{3)\2 2 z±^ + ^{ir(a* + p*)/2] sin[7r(a + (3)/2] 



sin[vr(/5+ (a + a*)/2)] 



x 



E(-i) 



m=0 



n + a 
m 



n + 3 
n — m 



x 



jr (-1)-' (^_^!_2m + 2m' 



m'=0 



n + a* 
m' 



n + 3 



n — m 
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X 



r(2±2 + n - m + m')T(^ + n + m - m') 
T(z±f + [3 + 2n) 



(17) 



This expression can also be evaluated directly relatively easily for the first few 
values of n. It also turns out that an a — a* can be factored out from the 
integral (17). 

From Eq. (2) we expect 



and we indeed get this expression by direct calculation from (15), (16) and 
(17). This expression vanishes (as it should) for real values of a, i.e. for un- 
broken PT symmetry. When the symmetry is broken [a is imaginary), then 
the second factor in (18) comes from the ratio of the gamma functions in the 
denominator of (17) and (16). 



6 Summary 

We have studied the Scarf II potential in its general form, which contains both 
its Hermitian and its PT symmetric version, depending on the choice of the 
potential parameters a and (3. Our work was motivated by the fact that this 
potential is a perfect laboratory to test the peculiarities of PT symmetric 
quantum mechanics, and in particular, the mechanism of the breakdown of 
PT symmetry. This is because it has no singularities on the x axis. 

The PT symmetric version of the Scarf II potential is obtained when (3 is 
chosen to be real, and a real or imaginary, depending on whether the PT 
symmetry is broken or unbroken: in the former case the bound-state energy 
spectrum is completely real, while in the latter case it is fully complex. The 
reality of (3 is not an essential restriction, since the roles of a and (3 can be 
interchanged in a trivial way. The potential is not sensitive to the sign of a, 
and this is reflected by the fact that both in the unbroken and broken PT 
symmetry case there are two sets of normalizable (bound) states with the 
same value of the principal quantum number n, and they are distinguished 
by the sign of a, i.e. the quasi-parity. The corresponding energy eigenvalues 
are different: they are both real for unbroken PT symmetry, and they form 
a complex conjugate pair when the symmetry is broken. 

We have devised a way by which the modified inner product of the wave- 



functions -?l"' /3 ' ±a ' ) = (V'n I -P|W ±a/3 ' ) ) can be evaluated explicitly. We have 



shown that these states are orthogonal in n and in the quasi-parity quantum 



Im(£(«./»)) 
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(a - a*)(a + a* + 2(3 + An + 2) , 



(18) 
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numbers, except when a is imaginary, i.e. the PT symmetry is broken, and 
the inner product of two states with the same n but different quasi-parities 
is considered. We derived a closed expression for the diagonal inner products, 
i.e. the pseudo-norm, and we found that its sign is indefinite, as expected from 
more general considerations. However, the sign of the pseudo-norm depends 
on a and j3 in a rather complicated way, so no compact rule can be formulated 
for its sign. 

Using the same techniques we also determined the standard inner product of 
the Hermitian Scarf II wavefunctions. In this case a* = (3 holds, and there is 
only one set of bound-state wavefunctions; the other one represents resonances: 
both the orthogonality of the states, the normalization constants for the Sacrf 
II potential have been derived for the first time. 

We also analysed the PT symmetric Scarf II potential as an ordinary complex 
potential. For this, we evaluated two types of diagonal matrix elements (using 
the standard inner product) of the PT symmetric wavefunctions: the matrix 
element of the imaginary component of the potential and the norm. The ratio 
of these gives a closed expression for the imaginaly component of the energy, 
which contains the potential parameters and thus can be used to study the 
breakdown of PT symmetry as the value of a is tuned from real values to 
imaginary ones through 0. In the case of other solvable PT symmetric poten- 
tials this transition requires crossing a singular point of the potential, which 
obscures important aspects of PT symmetry breaking. 



Acknowledgment 

This work was supported by the OTKA grant No. T031945 (Hungary), INFN 
and ENEA (Italy). 



Appendix 

Here we derive a general formula for the integrals which includes as special 
cases both the Hermitian and the PT symmetric inner product of the wave- 
functions. 

First, let's define the functions 

F^\x) = (1 - isinh(x))t+i(l + ismh(x))^Pj l a ' f3 \ismh(x)) , (19) 
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which are the bound-state Scarf II wavefunctions (3) without the normaliza- 
tion constant. The Jacobi polynomials appearing in (19) take the form [23] 



1 n 

P^(isinh(*)) = -£ 



m=0 




y n — m j 



x _ i S inh(x)) n - m (l + i sinh(:r)) m . (20) 

Also due to the properties of the Jacobi polynomials [23] the (19) functions 
transform under spatial reflection (i.e. the P operation) in the following way: 



F^\-x) = (1 + isinh(x))£+*(l - isinh(:r))i + 3P^)(-isinh(:r)) 

= (-l) n (l -isinh(a;))^ + 3(l + ismh(x))^Pjf' a) (ismh(x)) . 
= (-l) n F^ a \x) (21) 

Now define the integral 

oo 

Q { :{^ 5) = j F^\x)[F^ 5 \x)Ydx . (22) 



This integral contains the formulae necessary to evaluate the standard (Her- 
mitian) inner product of the ordinary (real) Scarf II bound-state eigenfunc- 
tions (with 7 = a and 8 — j3) and the PT symmetric inner product of the 
eigenf unctions of the PT symmetric Scarf potential with equal and different 
quasi-parity (with 7 = /3, 5 = a and —a, respectively). Note that in the latter 
case the roles of 7 and 5 are interchanged, due to (21). The above formula 
can also be used, of course, to evaluate integrals appearing in the Hermitian 
inner product of the eigenfunctions of the PT symmetric Scarf II potential, 
however, we shall find that in this case it cannot be reduced to simpler forms. 

First we need to evaluate integrals of the kind 



00 

J (l-isinhx) p (l + isinh:r) 9 da; , (23) 



which then appear in sums in (22). The integral can then be evaluated in a 
multistep procedure [20,27]: 



00 

A { ^ q) = j cosh p+9 xexp[(g-p)tanh~ 1 (isinha;)] (24) 
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J cosh p+l3 xexp[i(g — p) tan _1 (sinh:r)] 

-oo 
tt/2 

J cos~^ +<?+1 ^ y exp[i(g — p)y] 



-it/2 



F(-p - q) 



r(-p + l/2)r(-g + l/2) 



(25) 

(26) 
(27) 



In (24) we applied the relation tanh -1 z — \ In which is also used some- 
times to write the Scarf II wavefunctions in an alternative form [21,9,20], and 
in (26) a change of variables was performed: tany = sinhx. This method was 
used in Ref. [20] to evaluate the pseudo-norm of the ground-state wavefunc- 
tions for unbroken PT symmetry. In what follows we shall also find it useful 
to evaluate another integral, which follows directly from (27): 

oo 

A?' q) = J sinha;(l-isinha:) p (l + isinha;) 9 dx = i ^^~^ i 4 P '' ?) • (28) 



In the next step we have to evaluate sums of integrals of the type (23) appear- 
ing in (22): 



Q{aAi,8) = ( 1 yn+i 2 a+fi +f+ s *+2 s ' m l 7r (® + S*)/2] sm[n(P + 7*)] 




sin[7r(a + /3 + 7* + 5*)/2)] 
( 



m=0 



n + f3 \ ' 

E (-ir 

n — ml m'=o 



V 



m 



1 + 5* 
I — ml 



x 



V{^f +n -m + m' + l)T(^f- + l + m-m' + l) 
f( i±B±g + n + i + 2 ) 



(29) 



Without the loss of generality we can assume that n < I. With some rearrange- 
ment of the binomial coefficient and the gamma functions we can rewrite the 
last sum (over ml) into 



r(/ + 7 * + i)r(/ + £*) 
r( a+/3+7 * +5 * +n + l + l)n z 



E 



m 



+ n - m + ml + 1) Y(^f- + I + m - ml + 1) 



x ■ 



T(S* +m'+l) T(7* + l-m' + l) 

In what follows we assume that the relations 

5* = a 7* = /3 



(30) 



(31) 
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hold, in order to bring (30) to a simpler form. This assumption turns out 
to be valid for the most important cases, as it is described in the appropri- 
ate sections. Under the conditions (31) the two terms containing the ratios 
of gamma functions in (30) can be written as a finite power series of m'\ 
S"=o c j{ m 'y ■ This clearly follows from the two ratios of the gamma func- 
tions: the first one contains w! up to the (n — m)'th power in the form 
(a + m' + 1) (a + m' + 2) . . . (a + m' + n — m) , while the second one contains it 
up to the m'th power as {j3 + I — w! + l){j3 + I — w! + 2) . . . {j3 + 1 — w! + m). 
Now observing Eqs. 4.2.2.3 and 4.2.2.4 of Ref. [27] we find that 



E (-ir' E <h(™y = e c 3 e (-ir' I (™') J 

m'=0 \ m ' J 3=0 j=0 m'=0 \ w! 

n 

= E^ii(-l) , i! (32) 

3=0 

Since we assumed that n < I holds, j < I is also valid, so this sum can be 
non-zero only for n — l. In this case its value is 

5 ln (-l) n n\c n = 5 ln (-l) n n\(-l) m = 5 ln n\(-l) n+m . (33) 



Here we used that c n , the coefficient of (m') n (the highest possible power of 
m') in the power series is (+l) ra ~ m (— l) m , as can be seen from the factorization 
of the ratio of the gamma functions above. Substituting this result into (30) 
and (29) and then using Eq. 4.2.5.19 of Ref. [27], we get 



Q { y^' s=a * ] =6 nl (-iy 



a + 13 + 2n + 1 sin[7r(a + f3)\ 
i 





(34) 



Finally, we note that for the n = I diagonal case (29) can be evaluated in an 
alternative way too, using the sum rule 




■1)' 



n 
m 



(35) 



The interesting feature of this result is that the right hand side is independent 
of a and b. This formula is missing from the standard compilation [27], and 
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can be proven by an induction in n, using also the properties of the binomial 
coefficients. 



References 

[I] C. M. Bender, S. Boettcher, Phys. Rev. Lett. 24 (1998) 5243; 
J. Phys. A 31 (1998) L273. 

[2] C. M. Bender, S. Boetcher, P. N. Meisinger, J. Math. Phys. 40 (1999) 2201. 
[3] F. M. Fernandez, R. Guardiola, J. Ros, M. Znojil, J. Phys. A 32 (1999) 3105; 

C. M. Bender, G. V. Dunne, P. N. Meisinger, Phys. Lett. A 252 (1999) 272; 

M. Znojil, J. Phys. A 32 (1999) 7419. 
[4] F. M. Fernandez, R. Guardiola, J. Ros, M. Znojil, J. Phys. A 31 (1998) 10105; 

E. Caliceti, S. Graffi, M. Maioli, Commun. Math. Phys. 75 (1980) 51; 
C. M. Bender, G. V. Dunne, J. Math. Phys. 40 (1999) 4616. 

[5] M. Znojil, Phys. Lett. A 259 (1999) 220; 

J. Phys. A 33 (2000) 4561; 

J. Phys. A 33 (2000) L61; 

J. Phys. A 34 (2001) 9585. 
[6] G. Levai, M. Znojil, J. Phys. A 33 (2000) 7165. 
[7] B. Bagchi, R. Roychoudhury, J. Phys. A 33 (2000) LI; 

F. Cannata, M. Ioffe, R. Roychoudhury, P. Roy, Phys. Lett. A 281 (2001) 305. 
[8] F. Cannata, G. Junker, J. Trost, Phys. Lett. A 246 (1998) 219; 

B. Bagchi, S. Mallik, C. Quesne, Int. J. Mod. Phys. A 17 (2002) 51. 
[9] Z. Ahmed, Phys. Lett. A 282 (2001) 343. 
[10] M. Znojil, G. Levai, Mod. Phys. Lett. A 16 (2001) 2273. 

[II] A. Khare, B. P. Mandal, Phys. Lett. A 272 (2000) 53. 

[12] C. M. Bender, F. Cooper, P. N. Meisinger, V. M. Savage, Phys. Lett. A 259 
(1999) 224. 



[13] M. Znojil, |quant-ph/0103054| v3, 30 Oct. 2001. 



[14] G. S. Japaridze, J. Phys. A 35 (2002) 1709. 



15 



[15] L. E. Gendenshtein, Zh. Eksp. Teor. Fiz. Pis. Red. 38 (1983) 299 (Eng. transl. 
JETP Lett. 38 (1983) 35). 

[16] Z. Ahmed, Phys. Lett. A 290 (2001) 19. 

[17] B. Bagchi, C. Quesne, Phys. Lett. A 273 (2000) 285. 

[18] G. Levai, F. Cannata and A. Ventura, J. Phys. A 34 (2001) 839. 

[19] G. Levai, F. Cannata and A. Ventura, submitted. 

[20] B. Bagchi, C. Quesne and M. Znojil, Mod. Phys. Lett. A 16 (2001) 2047. 
[21] J. W. Dabrowska, A. Khare, U. P. Sukhatme, J. Phys. A 21 (1988) L195. 
[22] G. Levai, M. Znojil, Mod. Phys. Lett. A 16 (2001) 1973. 

[23] M. Abramowitz, I. A. Stegun, Handbook of Mathematical Functions, New York, 
Dover, 1970. 

[24] G. Levai, J. Phys. A 22 (1989) 689. 

[25] D. Baye, G. Levai, J.-M. Sparenberg, Nucl. Phys. A 599 (1996) 435. 

[26] A. A. Andrianov, M. V. Ioffe, F. Cannata, J. P. Dedonder, Int. J. Mod. Phys. 
A 14 (1999) 2675. 

[27] A. P. Prudnikov, Yu. A. Brychkov, O. I. Marichev, Integrals and Series, vol 1, 
Gordon And Breach, New York 1986. 



16 



